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^ ■ Abstract 

^ : 

The Greisen-Zatsepin-Kuzmin cutoff (GZK cutoff) predicted at the Ultra High Energy Cosmic 
'q^I Ray (UHECR) spectrum has been observed by the HiRes and Auger experiments. The results put 

O ' severe constraints on the effect of Lorentz Invariance Violation (LIV) which has been introduced 

^ . to explain the absence of the GZK cutoff indicated in the AGASA data. Assuming homogeneous 

a . 

source distribution with a single power-law spectrum, we calculate the spectrum observed on the 

(N ■ 

J> ' Earth by taking photopion production, e+e pair production, and the adiabatic energy loss into 

CN . account. The effect of LIV is also taken into account in the calculation. By fitting the HiRes 

i monocular spectra and the Auger combined spectra, we show that the LIV parameter is constrained 

(N ■ 

I ' to ^ = — O.S^Qg X 10~^^ and 0.01q4 x 10"^^ respectively, which is very consistent with strict Lorentz 

oo ; 

■ Invariance up to the highest energy. 

> : 
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I. INTRODUCTION 



A number of extensions of the standard model suggest that Lorentz Invariance (LI) is only 
a low-energy approximation and it may be deformed at very high energies, e.g., approaching 
the Planck Scale ~ lO^^ eV , y]. In a simple form shown by Coleman and Glashow, 

Lorentz Invariance Violation (LIV) can be expressed as a modified energy-momentum re- 



lation E'^ 



m 



^p^, under the assumption that LI is violated perturbatively in the 



context of conventional quantum field theory [s!]. It can also be interpreted in terms of 
different maximal attainable velocities for different particles. Modified energy-momentum 
relation with LIV terms proportional to the cube of the momentum, or a higher power, 
are also considered Although LI has been confirmed at accelerators up to 2 TeV for 
protons[7|], 104.5 GeV for electrons and 300 GeV for photons[8], it is still possible to see 
LIV in astrophysical processes with much higher energy, especially in the Ultra High Energy 
Cosmic Rays (UHECRs) with energies above 10^^ eV |9(. 

Since the gyration radius of UHECRs is larger than the height of our Galaxy in the 
Galactic magnetic field, UHECRs are generally thought to be of extragalactic origin. Prop- 
agating through intergalactic space, the UHECRs will interact with the Cosmic Microwave 
Background (CMB) photons, which results in energy and flux depletion. In particular, the 
photomeson process will induce a suppression in the spectrum above (3 — 6) x 10^^ eV and 
lead to the well-known Greisen-Zatsepin-Kuzmin (GZK) cutoff fiol . The spectrum of 
UHECRs can be calculated theoretically by assuming the source distribution and injection 
energy spectrum. The energy loss processes for UHECRs propagating in the intergalactic 
space include photopion production and e~^e~ pair production when interactin g w ith CMB, 
as well as the adiabatic energy loss due to the expansion of the Universe 



Be- 



cause of the short mean-free path of photopion production, the spectrum of UHECRs above 
6 X 10^^ eV falls sharply and results in the GZK cutoff. 

However, measurements of the UHECR spectrum have led to great confusion in the last 
decade. The result of the Akeno-AGASA experiment clearly shows an extension of the 



spectrum beyond the GZK cutoff 
cutoff, LIV has been introduced [5 
lead to the removal of GZK cutoff 



15|. To account for the AGASA data beyond the GZK 



. Even a LIV parameter as small as ^ ~ 3 x 10 



-23 



may 



16| . However, the measurements by HiRes 17| and 



Yakutsk p^] seem to show the existence of the GZK cutoff. 
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Recently, the situation tends to be clear. Having accumulated data for years, the HiRes 
Collaboration confirms their previous result and observes the GZK cutoff with a 5a standard 
deviation 19|]. The Pierre Auger Collaboration gives results consistent with HiRes and 



20 



2lj. As 



rejects a single power-law spectrum above ~ lO^^eV at the 6a confidence level 
confirmation of the GZK cutoff, severe constraints can be placed on the effect of LIV. 

In this paper, we investigate how the LIV can be constrained according to the latest 
HiRes and Auger data. We give details of the method to calculate the UHECR spectrum 
with LI in Sec. II. Then LIV is introduced to our calculation of the spectrum of UHECRs in 
Sec. HI. The modified spectrum with different LIV parameters and the constraints on the 
LIV parameters are also shown in this section. Section IV gives conclusions and discussions. 



II. THE SPECTRUM OF ULTRA HIGH ENERGY PROTONS IN THE STAN- 
DARD MODEL 

We assume the composition of UHECRs is pure proton. When propagating in intergalac- 
tic space, the ultra-high energy (UHE) protons will experience energy losses through the 
adiabatic expansion of the Universe, e~^e~ pair production, and photopion production due 
to the interaction with CMB photons. Then the energy evolution equation for a proton is 

-^^ = f^fiE) + Pf'^- (E) + P:{E) , (1) 

where l3f{E), (3^ ~ (E) and (3^{E) are the proton energy loss rate due to the Universe 
expansion, pair production and pion production respectively. 
The energy loss rates of protons at z = are jl^ . [3] 

f3f{E) = Ho, 



th 27 
rp poo 



2vr2^2 



de' a(e')K ( e')e' <- In 



th 



exp 



27T 



(2) 



where Hq is today's Hubble expansion rate, E and e are the energies of the proton and 
the CMB photon in the laboratory system (LS), respectively, e' is the photon energy in the 
proton rest system, 7 is the Lorentz factor of the proton in the LS, cr(e') is the interaction 
cross section, K{e') is the average fraction of energy loss, i.e., the inelasticity in the LS, n{e) 
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is the differential number density of CMB photons and T 2.73 K is the temperature of 
CMB. e(.i^ in Eq.Q is the threshold energy of the photon in the proton rest system above 
which the e^e~ pair or pion production can occur. Thus, ej.^^^ = 2me(l + me/mp) = 1.022 
MeV and e'^ = m,r(l + ^■w/'^f^p) = 149 MeV for p7 — > pe'^e~ and p7 pn respectively. 

Because of the redshifts of CMB photons, the energy loss rate will be larger at redshift 
z. The energy loss rate f3z{E) can be derived as 14 1 



/3fiE) = Hiz), /3r^-'-(E) = (1 + zfPoiil + z)E], 



(3) 



where H(z) = iJo\/^Af(l + z)^ + is the Hubble parameter. In this work we use the 
following cosmological parameters if n = 71 km s~^ Mpc~\ Qm = 0.27 and = 0.73 



according to the recent observations 23|. 

Solving Eq. ([1]) with the boundary condition E{z = 0) = Eq, we can get the initial 
energy distribution of protons which will be observed with energy Eq at the Earth. We 
denote this function as Eg{EQ, z), which means the initial energy distribution as a function 
of redshift z and observational energy Eq. We employ two assumptions: (1) proton sources 
are distributed homogeneously in the Universe without the evolution effect; (2) the source 
spectrum is a power-law with index jg. Then the observational proton spectrum at the 
Earth can be written as [l^ 



J{Eo) 



/•2max(-Eo) 


dt 


/ dz 


Jo 


dz 



(1 



z^Ep 



dEg{Eo,z) 



dEo 



with 

dEg{Eo,z) 
dEo 



[l + z) exp 



1 ^ (1+--)' 



dbojE'] 
dE' 



(4) 



(5) 



E'=(l+z')Eg{Eo,z')_ 

where 6o(^) = -dE/dt = E/3o{E), \dt/dz\ = 1/[(1 + z)H{z)], and ^max is the redshift of 
protons with maximum energy -Emax that reach us with energy Eq. In this work we adopt 
Emax ~ 10^^ eV. Larger -Emax do not affect the results. Lq is the total luminosity of UHE 
protons and is determined by matching the calculated spectrum to the observational data. 
(1 + z)"^ indicates the evolution of primary UHECR sources. However, this term is still 
unclear at present. Among the possible sources of UHECRs, the galaxies and some types of 
active galactic nuclei show m ^ 2.6 redshift evolution in radio, optical and X-ray bands 2J], 
while for BL Lacs there is a strong "negative" evolution 25j . It can be proven that for several 
reasonable evolution regimes, including the no-evolution case, the UHECR spectra can be 



reproduced well with different primary spectra[l^. Therefore, we adopt the no-evolution 
case for the source luminosity(m = 0) in this work. 

Eq-dlD can be simply understood: the protons within the initial energy interval {Eg, Eg + 
dEg) at redshift z contribute to the detected energy interval {Eq, Eq + dEo); the sum of all 
redshifts gives the total flux. 
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FIG. 1: Left panel: the expected UHE proton ener gy s pectrum (magenta curve) compared with 
the observational data from the AGASA HiRes 



19|], and Auger 



20( 1 experiments. The source 



spectrum of the theoretical curve is adopted as 7^ = 2.67, and the flux is normalized to the data 
of HiRes above 10^^ eV. Right panel: the fitting distribution as a function of the parameter 
7g for HiRes data. The three horizontal lines show the uncertainty ranges of 7^ at la, 2a and 3a" 
confidence levels, respectively. 



In Fig. [T] we show the theoretical spectrum of the UHE protons with the source index 
7c, = 2.67, which corresponds to the best fitting results of HiRes monocular data 13]. Data 
from AGASA 15| and Auger 20| are also shown in Fig. [H The theoretical flux has been 
normalized to the data of HiRes above 10^^ eV. The cosmic rays below 10^^ eV are usually 
thought to be of Galactic origin 26|], and we ignore them in this work. The right panel of 
Fig{T] shows the fltting distribution as a function of the source spectrum index 7^ using 
the HiRes data. The best fltting 7^ is 2.67 , with a x^/d.o.f. = 36.1/31. The statistical la 
range of 7^ is 2.58 < 7^ < 2.73. 
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III. THE SPECTRUM OF UHE PROTONS WITH LIV 



In this section we discuss how the LIV can modify the UHE proton spectrum. We 
adopt the framework of LIV developed by Coleman and Glashow, that a small first order 
perturbation is added to the free particle Lagrangian . The modification of the Lagrangian 
is then translated into the change of the dispersion relation of free particles, 

= + p'^ + ^p'^ , (6) 

where is the perturbation term. |^| is a very small parameter (~ 10^^^) and may be 
different for various particle species. 

Since the dispersion relation is changed by LIV, the kinematics of the reaction pj — >■ NX 
with N the nuclei and X the mesons, will also be changed, and finally, the UHE proton 
spectrum is modified. For a detailed treatment of the effect of LIV on the kinematics of 
P7 collision, we will closely follow the work of Alfaro and Palma 27|. In this work we only 
consider the process Ntt, for simplicity. 

The LIV effect generally modifies the kinematical inelasticity K, which is defined as the 
ratio between the energy of pion, E.^, and the primary proton energy Ep, K = E-^/Ep. The 



equation of inelasticity under the modified kinematics is given as 



27| 



(1 - Ke)V^ = F + j3cose^F^-SM{Ke) , (7) 

with 

s + SN{Ke) - s^{Ke) . . 

^ = vi ' 

where Kq is the inelasticity as a function of the angle 6 between the proton momentum 
in the center of mass system (CMS) and the direction of the CMS relative to the LS, 
s = E^^^ — pI^^ = {E + eY — {p — ky is the square of the total rest energy in the CMS, and 
and sn are the CMS energies of the pion and recoil nuclei, defined as 

Sa = El-pl = ml + ^aPl. (9) 

Note that for ^ < Sa can be negative and the particle will have a spacelikejour-momentum. 
We require Sa > in this work to guarantee that the particle is timelike [27,]. According to 
the definition of inelasticity, we have 

SN = ml + ^n[{1 - Ke)Ep]\ 

= ml + i^{KeEp)\ (10) 



in which we replace the perturbation term in Eq.Q with ^i?^, for simphcity. 



Following Ref. 



27l | we assume |^^| ^ |,^Ar| ~ 0, and denote as then Eq.Q can be 



further simplified. Eq.([7]) is solved numerically to get Kq, and the average with respect to 
9 gives the total inelasticity K = ^ Kgd6. We can see from the above equations that the 
inelasticity is a function of Ep and s = {E + e)"^ — {p — kY = Sp + 2y^e', where Sp is 
the energy of the proton in its CMS and e' is the photon energy in this system. For LI with 
= 0, -f^ is independent of Ep and is only a function of e'. 
In Fig. [2] we show the inelasticity as a function of the proton energy Ep in the LS and 
the CMB photon energy e' in the proton rest system, for the standard model as well as the 
models modified by LIV. Two kinds of LIV scenarios with ^ = 1 x 10^^'^ and ^ = — 1 x 10^^^ 
are investigated. We show that the modification of LIV is the reduction of the inelasticity 
at high energies. This effect is understood as the LIV leading to a reduction of the allowed 



phase space for the interaction 



27 



281]. It is interesting to note that no matter if ^ is positive 



or negative, the effect is similar — the inelasticity is reduced. Therefore, if the LIV exists, we 
can expect that the spectrum of UHE protons will be less suppressed by the jp interaction. 

The modified spectra of UHE protons for several values of LIV parameters are shown in 
Fig. [3], together with the unmodified spectrum from the standard model. We see that the 
GZK suppression effect becomes less significant for LIV cases. For very high energies or for 
large magnitudes of LIV parameters, the source spectra tend not to be distorted, i.e., the 
photopion production process '~fp —>■ Ntc does not play an important role any more. 

We employ a minimum fitting method to derive the implication of the HiRes and Auger 
data on the LIV parameter. A scan in the (^, 7^) plane is taken to calculate the UHECR 
spectra and the corresponding x^- Minimizing the distribution, we give a combined fit to 
get the source spectrum index 7^ and the LIV parameter ^ simultaneously. The confidence 
regions for la, 2a and 3a confidence levels in the ^ — 7g plane are shown in Fig. H] for the 
HiRes monocular spectra and the Auger combined spectrum, respectively. The best fitting 
results and la uncertainties of the parameters are 7^ = 2.671o;q2, ^ = ~^-8llo^ x 10^^^ for 
HiRes and 7c, = 2.57;^o;o2, ^ = ^-O-Ia ^ 10"^^ for Auger respectively. It is shown that the 
standard model with ^ = is very consistent with the present data. Compared to previous 
work, here we employ more strictly statistical analysis and include the negative part of ^. 

Considering an evolution factor (1 + z)^'^, we find that the best fitting results are 7^ = 
2.55tom, ^ = -O-ltli X lO'^s for HiRes data, and 7^ = 2.37;||;|]^, ^ = O.Oto.l x lO'^^ for 
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FIG. 2: The inelasticity as a function of the proton energy Ep in the LS and the CMB photon 
energy e' in the proton rest system, for the standard model (upper panel), LIV modified models 
with ^ = 1 X 10^^'^ (middle panel), and ^ = — 1 x 10~^^ (lower panel), respectively. 
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FIG. 3: The UHE proton spectra modified by LIV compared with the one predicted by the standard 
model. Left panel: > 0; right panel: ^ < 0. The observational data are from HiRes [iJ 
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FIG. 4: The la, 2a and 3a confidence regions (from inside to outside) on the £, ~ 7g parameter 



plane using HiRes (left panel [19:]) and Auger (right panel 
shows the best fitting values. 



3) 



) data above 10 eV. The red cross 



Auger surface data respectively. It is shown that the source spectrum differs a bit from the 
case with no source evolution. As for the LIV parameter, the results are very consistent 
within the statistical errors, and the conclusion is almost unchanged. 
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IV. CONCLUSION AND DISCUSSION 



The recent observations of the UHECR spectrum by HiRes [l9j and Auger [20|] show the 
existence of GZK suppression. In this work we use these observational data to test the LIV 
model and set constraint on the LIV parameter. The composition of UHECRs is assumed 
to be pure proton. We then solve the propagation equation of UHE protons in an expanding 
universe after incorporating the LIV effect in the energy loss rate of protons. A minimum 

fit to the observational data above 10^^ eV is adopted to derive the source spectrum of 
UHE protons and the LIV parameter. We find that the current data can limit the LIV 
parameter for pions to the level ~ 3 x 10"^'^. The standard model with the GZK cutoff is 
very consistent with the observational data. 

Only the LIV on the photopion production process is considered in this work. The pair 
production occurs at lower energy and the LIV effect might be less significant. To incorporate 
the LIV effect into the e^e~ production process, the analysis will be more complicated. In 
addition, we need to restrict the treatment to the case .^tt ^ ^at ~ 0. This has been shown 
to be due to the weakness of the presentation of the theory. We are unable to determine 
whether the perturbation term comes from the initial proton or the final state proton, which 
have different energy in the LS I27I. 

In this work the composition of UHECRs is assumed to be pure proton. However, the 



composition of UHECRs is poorly known from the experimental 



vations of HiRes show that the UHECRs are proton dominant 29|], while the results from 



point of view. The obser- 



the Auger experiment indicate a mediate mass composition [SOj. The determination of the 
UHECR composition depends on the interaction model and has a relatively large uncer- 
tainty at present. If the UHECRs are heavy nuclei dominant, the main process during the 
propagation in the CMB photon field is photo-disintegration, which will lead to the change 
of the conaposition and energy of primary cosmic rays and, accordingly, form a GZK-like 



spectrum [3l|, l32|, l33 |. 



We can make a rough estimate of the LIV effect in such a case, by assuming the primary 
composition of UHECRs to be iron and considering the process F e+''){C M B) — Mn+p. 
The analysis is greatly simplified by assuming ^Mn ~ ~ 0, while the physics is not 
changed as only the difference of the ^'s in the initial and final states is relevant for the 
kinematics j^. On the one hand, the interaction of photo-disintegration is possible only if 
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m|,g + AeEpe + ^FeEp^ > [niMn + mpY, resulting in ^pe > -4e2/[(mA/„ + m^)^ - mj,^] ^ 
—2 X 10^^^[e/eo]^ for CMB energy eo = 2.35 x lO^^eV. On the other hand, the spontaneous 
fragmentation ^^Fe — Mn + p should be forbidden for iron with energy lower than ~ 
3 X lO^^eV since cosmic rays with such energies have been detected. This condition requires 
■mp^ + C,FeEp^ < {mMn + fnpY, which gives ^pe < 1.2 x 10"^^ for Ep^ = 3 x 10^°eV. Therefore, 
we get —2 X 10"^^ < < 1-2 x 10~^^. It should be noted that the above estimates are quite 
rough. It will be more complicated to use the UHECR spectrum to study the LIV effects for 
heavy nuclei because there will be a chain of nuclei species taking effect in the interactions. 
Further detailed analysis is needed to probe the LIV if the UHECRs are proved to be heavy 
nuclei in future experiments. 
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